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We predict the polarization of A and A hyperons in peripheral heavy ion collisions at ultrarela- 
tivistic energy, based on the assumption of local thermodynamical equilibrium at freeze-out. The 
polarization vector is proportional to the curl of the inverse temperature field and its length, of the 
order of few percents, is maximal for particle with moderately high momentum lying on the reaction 
plane. A selective measurement of these particles could make A polarization detectable at the LHC. 
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In peripheral high energy heavy ion cohisions the sys- 
tem has a large angular momentum [1^. It has been re- 
cently shown in hydrodynamical computation that this 
leads to a large shear and vorticity [? . When the Quark- 
Gluon Plasma (QGP) is formed with low viscosity [3 , 
interesting new phenomena may occur like rotation [4 , 
or even turbulence, in the form of a starting Kelvin- 
Helmholtz instability (KHI) [S] [6], or other turbulent 
phenomena [7]. Furthermore, the large angular momen- 
tum may manifest itself in the polarization of secondary 
produced particles [11 [HI [9] . Recently, a formula for the 
polarization of weakly interacting particles with spin 1/2 
at local thermodynamical equilibrium has been found in 
Ref. [To] based on the extension of the Cooper-Frye for- 
mula to particles with spin. Provided that spin degrees 
of freedom equilibrate locally, the polarization turns out 
to be proportional to the vorticity of the inverse tem- 
perature four-vector field and can thus be predicted in 
a full hydrodynamical calculation of the collision process 
ended by the Cooper-Frye freeze-out prescription. 

Early measurements of the A hyperon polarization [TT] , 
averaged over a significantly large centrality range, in- 
dicated relatively small values, with an upper bound 
\Pj^ A I ^ 0-02 averaging over all azimuthal angles of A 
momentum. In this paper, we present a quantitative 
prediction of the A, A polarization, within a specific hy- 
drodynamical calculation, at different centralities at the 
LHC energy and its momentum dependence. Although 
the resulting polarization is of the order of 1% on aver- 
age (thus consistent with experimental RHIC bounds) , it 
turns out to be the largest (around 6-8%) for hyperons 
with moderately high momentum lying in the reaction 
plane. A selective measurement of As with few GeV mo- 
mentum into the reaction plane at the LHC could thus 
be able to show a finite polarization value, demonstrat- 
ing that also spin degrees of freedom achieve local equi- 
librium and, in an indirect way, that vorticous fiow is 
generated in peripheral heavy ion collisions. 

The A polarization arising from this pure thermo- 
mechanical effect (spin degrees of freedom equilibration 
due to the equipartition principle) in principle competes 



with the polarization induced by the electromagnetic 
fields. However, at the freeze-out stage, the magnetic 
field produced by the moving spectators is estimated to 
be of the order of 2 lO^^T [T? so that the resulting polar- 
ization is of the order of iin9kB/T ^ 10~^, i.e. at least 
four orders of magnitude less than our predicted value. 
The polarization of A hyperons has been approached with 
different models (e.g. [8[[T3]). Recently, Ref. [Il4] has con- 
sidered the local polarization of fermions in the plasma 
phase induced by the chiral anomaly, thus far with an un- 
specified transfer mechanism to final hadrons. We stress 
that in our approach the polarization of the observable 
hadrons is a consequence of the paradigm of local thermo- 
dynamical equilibrium; to be effective, the chiral anomaly 
should induce a modification of the velocity and temper- 
ature fields at the freeze-out. 



The A polarization in the participant centre-of-mass 
frame, as a function of its momentum, reads (in units 
c = i^ = l) ilOj: 
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where P^{x) — {1/T{x))u^{x) is the inverse temperature 
four-vector field, uf is the Fermi- Jiittner distribution of 
the A, that is l/(e^(^)-^-^(^) + l), being ^{x) = fi{x)/T{x) 
and p is the A four- momentum. Because at the temper- 
atures typical of freeze-out A is quite dilute {ttia ^ T), 
the Pauli blocking factor, (1— nj?), can be neglected in 
Eq. ([l]). The polarization vector is then proportional 
to the antisymmetric part of the gradient of the inverse 
temperature field, henceforth defined as thermal vortic- 
ity: 
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The spatial part of the polarization vector ([T]) gives rise 



to three terms: 
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The last two terms on the right hand side, involving po- 
lar vectors, should vanish because of the overall parity 
invar iance (achieved combining symmetry by reflection 
with respect to the reaction plane of the two colliding 
nuclei and invariance by rotation of tt around the axis 
orthogonal to the reaction plane). On the other hand, 
the first term, involving the spatial average of the curl 
of the (5 field, which is an axial vector, is not ought to 
vanish; in fact it is a vector aligned with the total an- 
gular momentum direction, which is orthogonal to the 
reaction plane (see Fig.[T]). It should be pointed out that 
these formulae apply to primary particles emitted from 
a locally equilibrated source. Secondary As emitted from 
either strong or weak decays - most likely - will have a 
lower polarization inherited from their parent particles. 
In the simplest scenario of an isochronous {t =const.) 
freeze-out at a given stage of the fluid dynamical expan- 
sion, according to Cooper-Frye prescription dT^xP^ -^ 
dVe, e = p^ being the A's energy. In this case, the above 
formula simplifies to: 



U{p) = 



he JdVnp (Vx/3) 



8m. 



JdVnp 



(4) 



The A polarization is usually determined by measuring 
the angular distribution of the decay protons, which, in 
the A rest frame is given by: 
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where a = 0.647, XIq is the polarization vector and p* 
is the direction of the decay proton, both in the A's rest 
frame. The vector XIq can thus be obtained by Lorentz 
boosting to this frame the one in Eq. Q: 



no(p) = n(p) 



e{e + m) 



n{p) ■ p 



(5) 



where (e, p) is A four-momentum and m its mass. One 
can readily realize that ||no|| < ||n|| and equality is 
achieved only if either p = (non-relativistic limit) or 
when p • n = 0. In both cases one has Tlo{p) = n(p). 
The above finding implies that maximal proper polariza- 
tion of A with finite momentum is achieved when they 
are transversely polarized. Thus, if n is directed along 
the total angular momentum {—y direction in Fig.fTl), A's 
having maximal polarization are those with momentum 
in the reaction plane or those with vanishing polar angle 
6 (normally undetectable) and, in this case, their proper 
polarization vector is aligned with the total angular mo- 
mentum. 
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FIG. 1. (Color online) Sketch of a peripheral heavy ion 
collisions at high energy. The A polarization points essentially 
into the direction of the total angular momentum {—y) of 
the interaction region, orthogonal to the reaction plane. As 
with the largest polarization are emitted into the (xz) reaction 
plane. 

The Hydro dynamical calculation: The goal of the 
hydrodynamic calculation is to evaluate the thermal vor- 
ticity ([2]) at the freeze-out. In this work we calculate it 
by using the Particle in Cell (PIC) fluid dynamic model, 
which provides us with the spacetime development of the 
flow of the QGP liquid until the freeze-out. In compari- 
son with Ref. [2^, only the relativistic case is considered. 

For computational purposes, it is convenient to absorb 
the h constant into (3^ and redefine thermal vorticity as: 



^M^ = -(S^/3^-5^/3^) 



(6) 



where /3^ = h/3^. Thereby, w becomes dimensionless. 
Note that in the thermal vorticity definition there is no 
projection of the derivatives transverse to the flow (the 
operator V^ = d^ — u^Uyd^)^ unlike in the usual defini- 
tion of the vorticity of the four- velocity field. 

We present in Fig. [2] the thermal vorticity val- 
ues weighted with the energy density, i.e. Vtzx = 
WzxEik/ {Eik). The weighting with the energy density 
of the cell is described in deteil in Ref. ^. From Fig. [2] 
it can be seen that, at the last time step presented, in 
the reaction plane we have already an extended area oc- 
cupied by matter. In case of peripheral reactions the 
multiplicity is relatively small, hence fluctuations in the 
reaction plane are considerable. In the relativistic case 
the outer edges show larger vorticity and random fluctua- 
tions are still strong. The average vorticity is smaller for 
the smaller impact parameters and it has positive value 
in the center and negative value at the edges. 

It should be pointed out that while the standard ve- 
locity field vorticity rapidly decreases with expansion [2j , 
thermal vorticity decrease is much slower and at some 
peripheral points it even increases. This is due to the 
fact that the matter cools during the expansion, so the 
temperature in the denominator of p^ decreases com- 
pensating for the decrease of velocity field vorticity with 
time. 



One should also mention that, in our calculation, hy- 
dro dynamical evolution starts after a dynamical longitu- 
dinal expansion based on collective Yang-Mills dynamics. 
The initial longitudinal size of the system is about 2x4 
fm, so the hydro process starts ~ 4 fm/c after the inter- 
penetration of the two Lorentz contracted nuclei. Conse- 
quently the configuration in Fig. [2] follows the interpen- 
etration by about 7.5 fm/c, which is the conventionally 
estimated freeze-out time. 
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FIG. 2. (Color online) The weighted average thermal vor- 
ticity, Qzx{x,z) (of the inverse temperature four- vector field 
/3^, see text for definition) calculated for all [x — z] layers 
at t=3.56 fm/c. The collision energy is ^/sNN = 2.76 TeV, 
b = 0.5bmax for (a) and b — 0.7 bmax for (b). Note that the 
configuration (a) does not favour KHI while (b) does. The 
cell size is dx = dy — dz — 0.585/0.4375 fm, while the aver- 
age weighted vorticity is {Q^zx) = 0.033 and 0.078 for (a) and 
(b) respectively. 



Results and Conclusion: The above described hy- 
drodynamical calculation was performed for the conser- 
vative case (a) presented in Fig. |2l which represents an 
initial rotation without the enhancement due to KHI. To 
calculate average polarization of A hyperons, the ther- 
mal vorticity has been properly weighted with the Fermi- 
Jiittner distribution np^ according to Eq. Q. 

The polarization vector, just as the flow vorticity, pri- 
marily points in the direction of the total angular mo- 
mentum {—y in Fig. IT]). It depends on the A's momen- 



tum vector through the Fermi- Jiittner distribution np^p) 
(see Eq. [4|. It increases with p^, and it is also sensi- 
tive to the flow properties and asymmetries; at ip^y ^ ^ 
GeV/c it is larger, about 6%, while in the reaction plane 
at Pa; ^ 3 GeV/c the polarization about 5%. The proper 
polarization vector in the A rest frame, determining the 
decay products angular distribution therein, is related to 
the polarization vector in the collision frame by Eq. (Isl), 
which introduces a further dependence on the A's direc- 
tion, as has been mentioned. Note that Eq. ^ modifies 
the direction of Ilo{p) with respect to 11 (p), except in the 
case when either p is lying on the reaction plane (py = 0) 
or orthogonal to the reaction plane {px^z = 0)- 
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FIG. 3. (Color online) The dominant ^/-component and 
the modulus of the observable polarization, IIo (p) in the rest 
frame of the emitted A as a function of the A momentum in 
the transverse plane at the participant cm. (i.e. at pz =0). 

The numerical results for the magnitude of the proper 
polarization and its projection along the angular momen- 
tum axis y are shown in Fig. (Is]). It can be seen that both 
increase as a function of transverse momentum and that 
their maximum values are attained in the reaction plane 



{x direction). The average value of polarization is of the 
order of 1% (consistent with RHIC bound), yet there are 
regions in momentum space where polarization is signif- 
icantly larger and could be found with a selective mea- 
surement. The maximum up to p^ > 3 GeV/c is around 
5%, when A's multiplicity is down by about 2 orders of 
magnitude compared to its top value, so that measuring 
polarization requires a sufficiently large statistics. These 
results are significant even if our hydrodynamical model 
assumes the possible maximum initial angular momen- 
tum; other (still realistic) models may have 10-20% less, 
yet without strongly reducing the final polarization value. 
It should also be noted that the same hydro model shows 
the possible occurrence of the KHI, which enhance the 
effect by 10-20% [5 . 

It is important to stress that, in order to measure po- 
larization, it is crucial to determine the orientation of 
the reaction plane, that is of the total angular momen- 
tum, on an event by event basis. As the polarization 
vector is oriented along the total angular momentum, a 
misidentification of the orientation would average to zero 
the measured polarization. A precise determination of 
the direction of the reaction plane is not as crucial be- 
cause A polarization does not vary much in length and 
direction around it (being at a maximum, see Fig. [3|. 
In order to improve accuracy, the participant center of 
mass should also be determined, both in pseudorapidity 
and in the transverse plane. This is usually not easy due 
to the limited acceptance of the central 47r detectors, but 
can be done by using the zero degree calorimeters with 
adequate correction factors as shown in Ref. [15], for the 
longitudinal cm. The same can be done in the transverse 
direction too. 



In conclusion, we have predicted the polarization of A 
hyperons in relativistic heavy ion collisions at the LHC 
energy and its momentum dependence. The polarization 
value may reach sizeable and detectable values of some 
percent in specific momentum regions, while the average 
value is predicted to be of the order of 1%, in agreement 
with the experimental bound set at RHIC energies [llj. 
In our hydrodynamical model, vorticity increases with 
the total angular momentum of the participant system 
in peripheral reactions, thus the observation of polariza- 
tion is more likely at the highest energy, that is at the 
LHC. Our calculation did not include the polarization 
of secondary As from decays of resonances or Ss which, 
most likely, will tend to dilute the signal. 

The observation of a polarization arising from this 
thermo-mechanical effect of equipartition of angular mo- 
mentum and in agreement with the predicted kinematical 
features would be a striking confirmation of the achieve- 
ment of local thermodynamical equilibrium (for the spin 
degrees of freedom too) of the matter created in relativis- 
tic heavy ion collisions. It would also indicate that signif- 
icant angular momentum predicted in [4 may persist up 
to the freeze-out. We stress again that this polarization 
exceeds by far the one induced by the magnetic field cre- 
ated by the spectators protons in the interaction region 

ini. 
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